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Abstract. In order to understand the parameter space ^^ of monic and centered 
complex polynomial vector fields in C of degree d, decomposed by the combinatorial 
classes of the vector fields, it is interesting to know the number of loci in parameter 
space consisting of vector fields with the same combinatorial data (corresponding 

^f^ ' to topological classification with fixed separatrices at infinity). 

^P ■ This paper answers questions posed by Adam L. Epstein and Tan Lei about 

the total number of combinatorial classes and the number of combinatorial classes 
corresponding to loci of a specific (real) dimension q in parameter space, for fixed 
degree d; these numbers are denoted by Cd and Cd,q respectively. These results 

f^ I are extensions of a result by Douady, Estrada, and Sentenac, which shows that 

the number of combinatorial classes of the structurally stable complex polynomial 
vector fields in C of degree d is the Catalan number Cd-i. 

We show that enumerating the combinatorial classes is equivalent to a so-called 
bracketing problem. Then we analyze the generating functions and find closed-form 
Cu I expressions for Cd and Cd,q, and wc furthermore make an asymptotic analysis of 

these sequences for d tending to cxj. 

These results arc also applicable to special classes of Abelian differentials, 
quadratic differentials with double poles, and singular holomorphic foliations of 
the plane. 



1. Introduction 

The space S^ of monic and centered single- variable complex polynomial vector fields 
of degree d is parameterized by its d— 1 complex coefficients. The space Ed — C^~^ 
is decomposed into combinatorial classes where the vector fields within each class 
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have the same combinatorial data set (to be defined). Each of these combmatorial 
classes is a connected manifold with well-defined (real) dimension q, which is the 
dim.ension of the com.binatorial class as a subspace in S^. 

Counting the combinatorial classes will help us to better understand S^, and 
the techniques utilized may also prove valuable to those who are considering other 
enumerative problems in dynamical systems. 

In this paper, we study the properties of the sequences Cd (the total number 
of combinatorial classes in S^) and Cd^q (the number of combinatorial classes of 
dimension q in 5^). These will be referred to as the simplified problem and the 
complete problem respectively. We first show in Section [3] that these problems 
are equivalent to so-called bracketing problems, where one counts the number of 
pairings of parentheses satisfying desired conditions [J . Then for each problem, a 
recursion equation and implicit expressions for the algebraic generating functions 
are calculated (Section [3. 4p . closed- form expressions for Cd and Cd,q are calculated 
(Section 2]), and asymptotic questions are considered (Section [5]). In particular, 
wc compute the asymptotic growth of the sequence Cd and prove that the discrete 
probability distribution iTd,q ~ ^^^ converges to a normal distribution for d — >■ oo. 
Lastly, we consider in Section [5] some related enumeration problems one might be 
interested in solving. 

Acknowledgements. Thanks to Adam Lawrence Epstein and Tan Lei for the 
suggestion of the problems and to Christian Henriksen for a simplification of the 
method for finding the generating functions. The author would furthermore like 
to thank Bodil Branner, Carsten Lunde Petersen, Carsten Thomassen, and again 
those mentioned above for helpful discussions and comments. The author is finally 
indebted to Philippe Flajolet for clarifying comments about his and his co-authors' 
work on analytic combinatorics. 

This research was supported by the European Union research training network 
CODY (Conformal Structures and Dynamics). 

2. Background and Definitions 

We present now a summary of some necessary concepts and definitions, and we 
include a new description of the combinatorial data set; it is different from the one 
presented in [2] (it is in fact more similar to the definition of the combinatorial 
invariant in [6]), but it is much more useful for our purposes here. For further 
details, please see [2] and [6]. 

It can be shown that oo is a pole of order d — 2 for vector fields ^p <E S^. There 
are 2{d — 1) trajectories 7^ which meet at infinity with asymptotic angles 2fd-i') ' 
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Figure 1. The point at oo is a pole of order d — 2 for vector fields ^p £ H^. There are 2{d ~ 1) 
trajectories 7^ which meet at infinity with asymptotic angles „, 7^ , , £ 6 {0, 1, . . . ,2d — 3}. When 
the labelling index i is even, the trajectories are called incoming to 00, and when the index £ odd, 

they are called outgoing from 00. 




Figure 2. There are 2d — 2 accesses to 00 defined by the trajectories at infinity. An end eg is 

infinity with access between 7f_i and 7^. An odd end is an end e^ labelled by an odd index k, 

and an even end is an end ej labelled by an even index j . 



i G {0, 1, . . . , 2d— 3}. When the labeUing index £ is even, the trajectories are called 
incoming to 00, and when the index £ odd, they are called outgoing from 00 (see 
Figure [1]). 

There are 2d— 2 accesses to 00 defined by the trajectories at infinity. An end ei 
is infinity with access between ji-i and 7^ (see Figure [2]). An odd end is an end Ck 
labelled by an odd index fc, and an even end is an end ej labelled by an even index 

J- 

Separatrices si are the maximal trajectories of ^p incoming to and outgoing from 
00. They are labelled by the 2{d—l) asymptotic angles 201^-1) ' ^ ^ i*^' l; ■ • ■ ' 2d— 3}. 
A separatrix s^ is called landing if si\ si = C, where C is an equilibrium point for 
^p. A separatrix Si = Skj is called homoclinic if Skj \ Skj = (see Figure |3l). A 
separatrix for a vector field ^p € S^ can only be either homoclinic or landing. A 
homoclinic separatrix Skj is labelled by the one odd index k and the one even index 
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Figure 3. Separatrices si are the maximal trajectories of ^p incoming to and outgoing from 
CO. They are labelled by the 2{d - 1) asymptotic angles ^^d-i) ' ^ ^ {0, 1, . . . , 2d - 3}. The 
separatrix structure can be equivalently represented in a separatrix disk model, by labelling the 
points cxp( 2d- 2 ^' £ = 0, . . . , 2d— 3 on §^ by sg and joining the points in the same equivalence class. 
In the figure, there are two homoclinic separatrices S3_o and sy^gi and four landing separatrices si, 
S2, S4, and S5. A few extra trajectories (in black) have been drawn in to better see the structure 
of the vector field. In this disk model, there are two sepal cells, one auj cell, and one center cell. 



j corresponding to its two asymptotic directions at infinity. The Separatrix graph: 

2d-3 

Tp = [J Si completely determines the topological structure of the trajectories of 

a vector field (see, for instance, |11| . [1]). 

In [2], the separatrix structure is encoded via an equivalence relation and a 
marked subset H on Z/(2(i — 2), where, in short, H tells which separatrices are 
homoclinic, and for the remaining separatrices (which are landing), the equivalence 
relation essentially tells which sg land at the same equilibrium point (the formal 
definition will be given in Subsection 12. 2p . 

This separatrix structure can be equivalently represented in a separatrix disk 
model, by labelling the points exp(^j3^), £ — 0, . . . ,2d— 3 onS^ by s^ and joining 
the points in the same equivalence class (see Figure [3]). 

Conversely, if we define an abstract combinatorial data set, (~, H), which consists 
of a non-crossing equivalence relation '^ and a marked subset H on Z/(2(i — 2) 
which satisfy certain properties (to be discussed), then it is known that there exists 
a monic and centered complex polynomial vector field of degree d whose separatrix 
structure matches the data set [2]. 

2.1. Zones The connected components Z oi C\ Tp arc called zones. There are 
three types of zones for vector fields in S^, and the types of zones are determined 
by the types of their boundaries: 

• A center zone Z contains an equilibrium point, which is a center, in its 
interior. Its boundary consists of one or several homoclinic separatrices and 
the point at infinity. If a center zone is on the left of n homoclinic separatrices 
Ski,ji , • ■ • , Sk„,j„ on the boundary dZ, then the center zone has n odd ends 
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Figure 4. Pictured are the trajectories of a vector field with four center zones; one odd center 
zone (shaded) with homochnic separatrices s^fi, si^2t and 83^4 and ends ei, 63, and 65 on the 
boundary; and three even center zones, each with one homoclinic separatrix and one end on the 

boundary. 



Ck-i , • • • , ek„ at infinity on dZ and the zone is called either a counter-clockwise 
center zone or an odd center zone. If a center zone is on the right of n 
homoclinic separatrices s^i ji , • ■ • , Sk„,j„ on the boundary dZ, then the center 
zone has n even ends e^^ , . . . , Cj^ at infinity on dZ and the zone is called either 
a clockwise center zone or an even center zone (see Figure [4]). 

• A sepal zone Z has exactly one equilibrium point on the boundary, which 
is both the a-limit point and w-limit point for all trajectories in Z (i.e. 
Ca = Cw)- This equilibrium point is necessarily a multiple equilibrium 
point. The boundary dZ contains exactly one incoming and one outgoing 
landing separatrix, and possibly one or several homoclinic separatrices, and 
the point at infinity. If a sepal zone is to the left of n homoclinic separatrices 
Sk-i,ji J • ■ • J Sknjn <3ii its boundary, then it has n + 1 odd ends on the boundary: 
Cfej^, . . . , efe„ and e^.+i for some corresponding j^ , depending on how one orders 
the separatrices. In this case, it is called an odd sepal zone. Similarly, if 
a sepal zone is on the right of n homoclinic separatrices Sk-^j-^, . . . , Sk„j„ 
on its boundary, then it has n+1 even ends on the boundary, e^j , . . . , Cj^ 
and Cki+i for some corresponding fc^, again depending on the ordering of the 
separatrices. In this case, it is called an even sepal zone (see Figure [5|). 

• An au-zone Z has two equilibrium points on the boundary, (a 7^ Ci^, the 
a-limit point and w-limit point for all trajectories in Z. The boundary dZ 
contains one or two incoming landing separatrices and one or two outgoing 
landing separatrices, possibly one or several homoclinic separatrices, and 
the point at infinity. If an aw-zonc is both on the left of ni homoclinic 
separatrices Sk-^j-^ 1 ■ ■ ■ i Sk„ ,j„ and on the right of n2 homoclinic separatrices 
Ski.jn • ■ • , s/c„ ,jn on the boundary, then the aw-zone has rii + 1 odd ends 
(efcj^ , . . . , efc^ and ej.-^-i for some corresponding ji) and 712 + 1 even ends 
(cjj^ , . . . , Cj^ and Ck^+i for some corresponding ki) on the boundary (see 
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Figure 5. Pictured arc the trajectories of a vector field with an even sepal zone (shaded). On 
the boundary of the sepal zone is the double equilibrium point which is both the a and ui limit 
point of the trajectories; one incoming landing separatrix 34 and one outgoing landing separatrix 
S5; three hoinoclinic scparatrices sifi, 83^2, and st^q; and four ends at infinity eo, £2, 64, and 
66- There is an odd sepal zone (not shaded) which shares the equilibrium point and the landing 
separatrices with the shaded sepal zone, but it has no homoclinic scparatrices and only one odd 

end 65 on the boundary. 




Figure 6. Pictured are the trajectories of a vector field with an aa;-zone (shaded). On the 
boundary of the zone are two equilibrium points: one which is the a-limit point of the trajectories 
and the other is the oj-limit point of the trajectories. Also on the boundary are one incoming 
landing separatrix sq and one outgoing landing separatrix S3. The zone is to the left of the 
homoclinic separatrix si_2 and on the right of the two homoclinic scparatrices 55,4, and S7,6- 
Finally, the boundary contains two odd ends ei and £3 and three even ends eo, 64, and eg at 

infinity. 



Figure ^ . 

Remark. It will be important to note that for a sepal zone, there is exactly one end 
whose index is not equal to an index of a homoclinic separatrix (in the notation 
above, eji_(.i for an odd sepal zone and ek^+i for an even sepal zone). Similarly, 
for an aw-zone, there are exactly one odd end and one even end, neither of whose 
indices coincide with any index of a homoclinic separatrix (in the notation above 
the odd and even ends are e^.+i and Ck^+i respectively). 
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In the separatrix disk model, the connected components of B \ IJ s^ are called 

e 
cells. Cells are the abstract analog of zones, and are named accordingly (i.e. there 

are center cells, sepal cells, and atu cells). The types of cells are determined by the 

types of their boundaries (see Figure [3]). 

2.2. Combinatorial Data Set Now that we have defined the types of zones and 
corresponding cells, we can state the definition of the combinatorial data set as 
presented in [2]. 

Definition. A combinatorial data set {^,H) of degree d > 2 consists of an 
equivalence relation ~ on Z/ {2d — 2) and a marked subset H C "Z/ {2d — 2) 
satisfying: 

1) '-^ is non-crossing. 

2) If £' y^ £", then £' - £" and i' e H ^ e" e H and f and £" have different 
parity. 

3) Every cell in the disk-model realization of (~, H) is one of the five types: an 
aw-cell, an odd or even sepal-cell, or an odd or even center-cell characterized 
as above. 

2.3. Transversals There is an equivalent way to encode the combinatorial 
structure of a vector field, which is much more useful for our purposes here. We 
define in this section the important structures needed to understand this equivalent 
definition of a combinatorial data set. 

In any simply connected domain avoiding zeros of P, the differential -^|y has an 
antiderivative, unique up to addition by a constant 

Note that 

cj,^{^p)^^'{z)P{z)±^±. (1) 

The coordinates w — 4>{z) are, for this reason, called rectifying coordinates. We will 
call the images of zones under rectifying coordinates rectified zones. The rectified 
zones are of the following types: 

• The image of an aw-zonc under (f> is a, horizontal strip (see Figure [7]). 

• The image of an odd sepal zone under is an upper half plane, and the image 
of an even sepal zone is a lower half plane (see Figure |8]). 

• The image of a center zone (minus a curve contained in the zone which joins 
the center ( and cxd) under cp is a. vertical half strip. It is an upper vertical 
half strip for an odd center zone and a lower vertical half strip for an even 
center zone (see Figure [S]). 
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Figure 7. The image of an aixi-zonc under </> is a horizontal strip. 
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Figure 8. The image of an odd sepal zone under <f> is an upper half plane, and the image of an 
even sepal zone is a lower half plane. In this figure, there is an even sepal zone mapped to a lower 

half plane. 



Via the rectifying coordinates, it is evident that there are a number of closed 
geodesies in C \ {equihbrium pts} in the metric with length element ipAi through 
00. Among these are the h homoclinic separatrices, and there are s distinguished 
transversals (defined below). 

Definition. The distinguished transversal Tkj is the geodesic in the metric rpr^i 
joining the ends e^ and Cj, avoiding the separatrices and equilibrium points, where 
Cj is the left-most end on the upper boundary and e^ is the right-most end on the 
lower boundary of the strip which is the image of the aw-zone the transversal is 
contained in (see Figure [TOl and [TT]) . 

Note that the way in which the distinguished transversal is chosen, the indices 
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Figure 9. The image of a center zone (minus a curve contained in the zone which joins the center 

f and 00) under </) is a vertical half strip. It is an upper vertical half strip for an odd center zone 

and a lower vertical half strip for an even center zone. In this case, there is an odd center zone 

mapped to an upper vertical half strip. 
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Figure 10. There may be several transversal which avoid the equilibrium points and separatrices 
(the dashed curves), but there is exactly one distinguished transversal for each atj-zone (in this 

case, T-ifi). 



of the ends it joins are exactly those ends whose indices will never coincide with 
the indices of any homoclinic separatrices. 

2.4. Alternative Combinatorial Description We will prove that one can 
equivalently describe the combinatorics by the union of homoclinic separatrices 
Sk,j and distinguished transversals Tj-j. We state the definition below, and prove 
that it is equivalent to the one presented before (from [2]). 

Definition. A combinatorial data set (i, H, T) of degree d > 2 consists of an 
involution t : Z/(2(i — 2) — > 7^/ {2d — 2) and marked subsets H and T satisfying: 

• If l{£) ^ i, then £ and i{£) have different parity (if one is even, then the other 
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Figure 11. Each acj-zone is isomorphic to a strip. We define the distinguished transversal to be 
the geodesic in the metric |(i2:|/|P(2)| joining the ends ej; and ej (in this figure, 63 and eo) where 
Bj is the left-most end on the upper boundary of the strip and e^ is the right-most end on the 
lower boundary of the strip. Since the indices of the distinguished transversal are the same as the 
indices for the two landing separatrices on the upper left and lower right boundary of the strip, 
they can never coincide with the indices for a homoclinic separatrix. 



is odd). 

• The involution l is non-crossing. 

• If i{£) y^ £, then cither both £ and l{£) arc in H or both £ and l{£) are in T. 

This definition leads to a description of the combinatorics via non-crossing 
pairings of even and odd numbers. The pairing of indices in H represents pairings 
of indices k and j corresponding to homoclinic separatrices Sk.j, and the pairing 
of indices in T represents pairings of indices k and j corresponding to indices of 
ends efe and Cj that the distinguished transversal Tj^j joins. The indices which 
are not paired (i.e. £ such that i,{£) = £) correspond to the indices of ends ef on 
the boundary of sepal zones which do not coincide with indices of any homoclinic 
separatrices. Essentially, we want to use the numbers Z/(2(i — 2) to stand for 
indices of separatrices for homoclinics, and indices of ends otherwise. Since these 
were chosen in the way described above, this never causes a conflict. 

As one can draw a disk model with a separatrix graph given the first definition 
of combinatorial invariant, one can draw a transversal graph given the second 
definition. The transversal structure can be equivalently represented in a transversal 
disk models by labelling the points exp( „^?^^„ ), £ = 0, . . . , 2(i — 3 on S^ by S£ and 



the points between them exp( 



27ri(f-l/2) 
2d-2 



) by Bf, joining sj. and Sj corresponding 



to homoclinic separatrices Sfej. and joining Ck and Cj corresponding to each 
distinguished transversal T^j (see Figure fT2|) . That the involution l is non-crossing 
means that the transversal and homoclinic separatrices in the transversal disk model 
do not cross. 
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Figure 12. Pictured is the transversal disk model equivalent to the separatrix disk model pictured 

in Figure [3] There are two homoclinic separatrices S3,o and S7_6, and there is one distinguished 

transversal T5,4. Note that each connected component contains exactly one equilibrium point. 



Proposition 2.1. The definition above is equivalent to the definition of 
combinatorial data set presented in f^ (Definition \2.2\) 

Proof. We basically need to prove that separatrix disk raodels are in one-to-one 
correspondence with transversal disk models. 

If we start with a separatrix disk model, such that every cell is one of the five 
types, then we can construct a unique transversal disk model. Note first that the 
homoclinic separatrices stay the same. Then for each aw-zone, there are exactly 
one even end Cj and one odd end e^ whose indices do not coincide with any of the 
homoclinic indices on the boundary. Join these ends in the disk. Finally, remove 
all landing separatrices from the separatrix disk model, and what you are left with 
is a transversal disk model (see Figure [T^ . 

Now conversely given a valid transversal disk model (i.e. one satisfying the 
definition above), we again leave the homoclinic separatrices alone as before. The 
homoclinic separatrices and distinguished transversals decompose the transversal 
disk model into /i + s + 1 connected components. Call each of these components 
a transversal cell (as not to confuse them with connected components with respect 
to a separatrix disk model) . Let all Si on the boundary of each transversal cell be 
joined to the same point in the interior of the transversal cell. If there are no S( 
on the boundary of a transversal cell, then that component is a center cell (with 
respect to the separatrix disk model) . Now remove the curves corresponding to the 
distinguished transversals. This procedure gives a separatrix disk model, and now 
we just need to show that this separatrix disk model has the allowable types of cells. 
As already mentioned, if there are only homoclinic separatrices on the boundary 
of a separatrix cell, then that component is a center cell. If there is an open chain 
of n (where we allow n = 0) homoclinic separatrices: Sk-^j-^, ■ ■ ■ ,Sk„,j„ such that 
either ji^i = ki + 1 for z = 1, . . . , n — 1 and Sji_i and Sfc„+i are landing (see Figure 
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Figure 13. An open chain of three homoclinic separatrices sj.^ j^, s^^j^, and s^^j^ such that 
ji-\-i = fci + 1 for j = 1, 2 and Sj-^ — i and sjjg+i are landing. 




Figure 14. An open chain of three homoclinic separatrices s/s^^jj, Sfe2,J2 

fci.^! = j'i + 1 for i = 1, 2 and sj;j_i and Sj„+i are landing 



and Sfe3,J3 such that 



[T3D or if k 



i+l 



h + 



1 for i = 1, 



1 and Sfci-i and Sj„+i are landing (see 



Figure [H)) on the boundary of the same separatrix cch, then either the two landing 
separatrices land at the same point, or they do not. If they do land at the same 
point, then Sjj_i and Sfc„+i (or s^i-i and Sj„+i) are the two landing separatrices 
on the boundary of a sepal cell. If they land at different points, then they must be 
separated by a distinguished transversal Tk^k^+i for some k associated to another 
open homoclinic chain (or Tj^+i j for some j associated to another open homoclinic 
chain) on the boundary of the same separatrix cell. Then Sj-^-i must land at the 
same point as Sfc, since they arc on the boundary of the same connected component 
in the transversal disk model (similarly for Sfcj-i and Sj) (see Figure [T5|) . One can 
make a similar argument for the other open homoclinic chain involved, and then 
one must get an aw-cell. □ 



3. Setup of the Bracketing Problem 

We have seen that the combinatorial structure of a monic and centered complex 
polynomial vector field can be completely described by pairings of even integers with 
odd integers, and conversely that any such pairings as defined by the transversal 
combinatorial data set leads to the combinatorics for a monic and centered complex 
polynomial vector field. The goal is to utilize the alternative combinatorial 
description in order to convert the problem of counting the combinatorial classes 
into a so-called bracketing problem: a combinatorial problem involving pairings of 
parentheses placed in a string of elements (see [4]) in a valid way (to be defined). 
In short, each homoclinic separatrix Sk,j has exactly one even and one odd index 
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Figure 15. An open chain of three honioclinic separatrices Sfcijj, •5fe2,J2' 
ki^l = ji + 1 for j = 1, 2 and Sfej_i and Sj„+i are landing. If they land at different points, then 
they must be separated by a distinguished transversal Tjj+ij for some j associated to another 
open homoclinic chain on the boundary of the same separatrix cell. Then Sfej_i must land at 
the same point as Sj, since they are on the boundary of the same connected component in the 

transversal disk model. 



associated with it, and each distinguished transversal T^j has exactly one even and 
one odd index associated to it. Therefore, each combinatorial data set can be fully 
described by pairings of even and odd numbers, corresponding to these mutually 
disjoint pairs of indices, and hence they can be counted by a bracketing problem. 
In our case, we will have balanced parenthetical configurations placed in the string 
1 2 . . .2d — 3, representing the indices in Z/{2d — 2), and we will use round 
brackets (...) to pair an even an odd number corresponding to homoclinics and 
square brackets [. . . ] to pair an even and odd number corresponding to distinguished 
transversals. We now elaborate on these ideas. 

3.1. Structurally Stable Vector Fields For peadagogical reasons, we start by 
describing the case of counting structurally stable complex polynomial vector fields 
by a bracketing problem. 

Theorem 3.1 (Douady, Estrada, Sentenac) The number of combinatorial 
classes for the structurally stable vector fields in S^ is the Catalan number 

Proof For structurally stable vector fields, the combinatorics is completely 
describable by the d—1 transversals (since there are neither homoclinic separatrices 
nor sepal zones), and hence non-crossing pairings of each odd end e^ with an 
even end ej (this is, in fact, how j^ defines the combinatorial invariant for the 
structurally stable vector fields) . This pairing of ends can be represented by placing 
parentheses in the string 1 2 ... 2d— 3, where the elements paired by parentheses 
correspond to the labels of the ends wc want to pairp Note that this is equivalent 

f A more familiar description might be evident in the disk model, where this pairing can be seen 
as the number of non-crossing handshakes of 2{d — 1) people seated around a round table. 
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to the number of ways to make valid pairings of d — 1 pairs of parentheses, where 
the string 1 2 ... 2d — 3 is not actuahy necessary, since every end is paired with 
another. For example, if d — 1 = 3, then the configurations in the table below are 
equivalent: 



[01] [2 [34] 5] 


[][[]] 


[0 [12] 3] [45] 


[[]][] 


[0[1[23]4]5] 


[[[]]] 


[0 [12] [34] 5] 


[[][]] 


[01] [23] [45] 


[][][] 



You might recognize then that the number of combinatorial classes (with the 
labelling of the scparatriccs) for the structurally stable monic and centered complex 
polynomial vector fields of degree d is just the Catalan number Cd-i- The proof is 
classical (see for instance |5]), but we include it here for completeness. 

The goal is to write a recursion equation for C„, the n**^ Catalan number, which 
counts the number of configurations matching n pairs of parentheses. The first 
character of any balanced configuration is an open parenthesis " [" . Somewhere 
in the configuration is the matching "]" for the open one. In between that pair 
of parentheses is a balanced configuration, and to the right is another balanced 
configuration: 

[A]B, (3) 

where A and B arc balanced parenthetical configurations. Each A and B can have 
anywhere from to n — 1 pairs of parentheses, but together they must have exactly 
n — 1 pairs of parentheses. So if A has k pairs of parentheses, then B must have 
n — fc— 1 pairs of parentheses. Thus we count all configurations where A has pairs 
and B has n — 1 pairs, A has 1 pair and B has n — 2 pairs, etc. Add them up, and 
you get the total number of configurations with n balanced pairs of parentheses. 
Therefore, 

Cn = CoCn-l + CiCn-2 + ' ' ' + Cn-'lCl + C„_iCo, Cq = 1. (4) 

We will use this recursion equation and the generating function 



/(Z) = J2 Cn^ 



(5) 



71=0 



in order to derive a closed- form expression for C„. Using ^, one can compute 

[/(z)]' = CoCa + {CoCi + CiCo) z + (C0C2 + CiCi + C2C0) z^ + .... (6) 
Combining this with (j?]), one arrives at 



Therefore, 



[f{z)f = Ci+ C2Z + C32' + Ciz^ + 



f{z)=C^ + z[f{z)Y 



(7) 
(8) 



Prepared using etds.cis 



Enumerating combinatorial classes of vector fields 



15 




65 

[0 1][2[3 4]5] 



So S3 



Co Ei 




ei es 



So S3 



Co £4 




[0[12]3][4 5] 



[0[1[2 3]4]5] 




ei 63 



[0[12][3 4]5] 




[0 1][2 3] [4 5] 



Figure 16. The five different disk models for the structurally stable vector fields of degree d = 4. 

The involution on the ends is marked by the dashed curves. The representation of the involution 

in brackets is displayed below each figure. 



Applying the quadratic formula, 
and one must choose the branch 



2z 



1- Vl~4i 

2z 



(9) 



(10) 



since f{z) — > 1 as z — > 0, which is desired since f{z) = Cq = 1 (choosing the "+" 
branch gives f{z) — > oo as z — >■ 0). Expanding ^/1^^4z = (1 — 4z) using the 
binomial formula and making some arithmetic manipulations, one arrives at 



n=0 ^ ^ 

giving the desired result. 

An example is given in Figure [16] to demonstrate. 



(11) 



n 



3.2. Nan- Structurally Stable Vector Fields The convention, as stated before, is 
that round parentheses (• • • ) are used to mark pairings corresponding to homoclinic 
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[0 1][2 3 4 5] 



0[1 2]3[4 5] 



Figure 17. Disk models for two examples of vector fields of degree d = 4 having sepal zones 

and no homoclinic scparatriccs. The pairing of the ends is marked by the dashed curves. The 

representation of the combinatorics in brackets is displayed below each figure. 



separatrices, square parentheses [• • • ] are used to mark pairings corresponding to 
a distinguished transversal in each aw-zone, and elements that are not paired 
correspond to the sepal-zones. 

We first look at the case when there are sepal zones but no homoclinic 
separatrices. Therefore, the zones can only be cither aw-zoncs, having exactly 
one even and one odd end, or sepal zones, having exactly one even or one odd end. 
The string 1 2...2d—3 represents the labels of the 2d— 2 ends in this case. The 
aw-zones will then be denoted by square parentheses [• • • ] in this string, and the 
labels corresponding to the ends of sepal zones are not paired by a set of parentheses. 
We see now that the string 1 2 ... 2d — 3 is necessary since the unpaired elements 
will be placeholders corresponding to the sepal zones. Two examples are given in 
Figure [m 

What then happens if we introduce homoclinic separatrices? Then it is no longer 
enough to consider pairings of ends. The string 1 2 . . .2d — 3 now represents 
the indices which may belong to either the ends ei or separatrices Sf (as was in 
the involution definition of the combinatorial data set). For every homoclinic 
separatrix Sk.j, we pair the numbers k and j by round parentheses (■■■)■ For 
example, 1 . . . {j . . . k) . . . or 1 ... (A; ... j) ... . As usual, for every distinguished 
transversal Tkj , we pair the numbers k and j by square parentheses [• • • ] (see Figure 
[HD. 

We can see from this that given any combinatorial data set, there is a unique 
bracketing representation defined in this way. 

3.3. Valid Bracketings The bracketing must satisfy certain rules so that they are 
in accordance with what can happen for a vector field (and hence a combinatorial 
data set). Pairs of square and round parentheses placed in a string of elements 
1 2 ... 2(i — 3 is called a valid bracketing if: 

1. there are an equal number of right and left parentheses for each type, 

2. the number of left parentheses must be greater than or equal to the number 
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ei 6.1 



[0(1 2)3] [4 5] 




(0[1(2 3)4]5 



Figure 18. Disk models for three examples of vector fields of degree d = 4 having sepal zones 

or/and homoclinic separatrices. The pairing of the ends is marked by the dashed curves. The 

representation of the combinatorics in brackets is displayed below each figure. 



of right, reading from left to right Example: ())()( not valid, 

3. there must be at least one element between successive left (resp. right) 
parentheses Example: ((01)23) not valid, 

4. there must be an even number of elements in each pair of parentheses 
Example: (012)3 not valid, 

5. and the square brackets must not cross the round brackets Example: ([)] 
not valid 

It is clear that a valid bracketing as defined above leads to an admissible 
combinatorial data set. 



3.4. Recursion Equations and Generating Functions We compute the recursion 
equation for p„, the number of ways to place any number of round or square 
parentheses in a valid way in a string of n elements, by using the following algorithm 
to generate any possible valid parenthetical configuration in a string of n elements, 
which we will denote by s„. 



's"s„_i if string does not start with parenthesis 

'(s"s2a"s)"sh if string starts with round parentheses (12) 

" [s"s2a"s]"sh if string starts with square parentheses 



for non-negative integers a and b satisfying 2a + 6 + 2 = n, and the s representing 
numbers in {0, 1, . . . , 2d — 3}. A table of the construction of such strings from n = 1 
to 71 = 5 is given below. 
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n 


S Sn-1 


a 


b 


"(s"S2a"s)"Sfc 


"[s"S2a"s]"S6 


1 


S 










2 


ss 








(ss) 


[ss] 


3 


sss 

s{ss) 

s[ss] 





1 


{ss)s 


[ss]s 


4 


ssss 





2 


{ss)ss 


[ss]ss 




ss{ss) 






{ss){ss) 


[ss]{ss) 




ss[ss] 






{ss)[ss] 


[ss][ss] 




s{ss)s 


1 





(ssss) 


[ssss] 




s[ss]s 






{s{ss)s) 
{s[ss]s) 


[s{ss)s] 
[s[ss]s] 


5 


sssss 





3 


{ss)sss 


[ss]sss 




sss{ss) 






{ss)s{ss) 


[ss]s{ss) 




sss[ss] 






{ss)s[ss] 


[ss]s[ss] 




ss{ss)s 






{ss){ss)s 


[ss]{ss)s 




ss[ss]s 






{ss)[ss]s 


[ss][ss]s 




s{ss){ss) 


1 


1 


{ssss)s 


[ssss]s 




s[ss]{ss) 






{s(ss)s)s 


[s{ss)s]s 




s{ss)[ss] 






{s[ss]s)s 


[s[ss]s]s 




s[ss][ss] 












s(ssss) 












s[ssss] 












s(s(ss)s) 












s[s(ss)s] 












s(s[ss]s) 












s[s[ss]s] 











Then p„ satisfies the recursion equation 



Pn-l +2 Y. P2aPb 

\2a+b+2=n . 

1, 

0, 



n > 1 and a, 6 > 



n < 0. 



(13) 



We are interested in the generating function 



G(z)=^CdZ^ 



d=\ 



(14) 



where c^ = P2(d-i)5 since for degree d the string wc work with is . . . 2(i — 3. Note 
that 0(0) — 0. For n > 0, let g„ = -p^n and r„ = Vin-\ {tq = P-i = 0). Let also 



/(2) = XI 9"" 



(15) 



n^O 



Prepared using etds.cis 



Enumerating combinatorial classes of vector fields 19 

and 

.g(z) = 5^r„z". (16) 

n = l 

Note that zf — G. One can deduce for n > 

n-l 

qn = P2n = P271-1 + 2 ^ P2]P2n-2]-2 
3=0 
n-1 

= rn + 2^q.jqn-j^i (17) 

3=0 

and 

rn = P2n-1 = P2n-2 + 2 ^ P2jP2n-2j~3. 
n-1 

= (7„_i + 2 ^ q.jrn-j-1. (18) 

i=o 

From equations p7)) and (|18p . one can deduce 

/ = l + .9 + 2z/2 and g = z/ + 2z/5 (19) 

respectively, giving 

4zV'-4z/2 + (z + l)/-l = 0. (20) 

Substituting z/ = G, one gets that G satisfies the algebraic equation 

G3_G2 + ^i±^G-- = 0. (21) 

4 4 ^ ^ 

For the complete problem, namely giving the number of combinatorial classes of 
each dimension, we are interested in the generating function 

oo 
G{z,t)= ^ Cd,qzH'', Cd,q=P2{d-l),q, (22) 

d=l,g=0 

where q is the real dimension of the combinatorial class in parameter space. Note 
that 

oo 

G(z,l)= Y. Cd,qz'' = G{z). (23) 

d=l,q=0 

FoUlowing the same method as before, 

's"s„_i^q if string does not start with parenthesis 

-",,- ^^ '(S"s2a,gi"s)" 56,92 91+92 + 1=9 (24) 



"[s"S2a,9i"s]"Sfc^,, qi+q2 + 2 = q 



for 2a + 6 + 2 ~ n. The condidtions on qi and (72 conic from the fact that 
each homoclinic separatrix (and hence each round pair of parentheses) contributes 
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real dimension 1, and each aoj-zone (and hence each square pair of parentheses) 
contributes complex dimension 1 (real dimension 2). Manipulating the coefficients 
as in equations (|17l) and (IT51) . we get the equations 

f^l+g + z{tf+t^f) and g ^ zf + tzfg + t^zfg. (25) 

We then have that G = zf satisfies the algebraic equation 

[t + efG^ - 2(t + t2)G2 + (1 - z + z(f + e))G - z = 0. (26) 



4. Closed Forms 

The best one can hope to achieve with respect to an enumerative problem is to find 
a closed-form representation for that what is being enumerated. The main tool we 
will use is the Lagrangc-Biirmann Inversion Theorem. We will only use a simplified 
version and state it here for completeness. 

Theorem 4.1 (Lagrange-Burmann Inversion Theorem) Let 0(m) he a for- 
mal power series with (J)q ^ 0, and let Y(z") he the unique formal power series 
solution of the equation Y = z(f){Y). Then the coefficient ofY{z) of order n is 
given by 

ny(z) = iK-i]0«. (27) 

The notation [z"] • • • means the coefficent of z" in the power series expansion of 
the expression that follows. 

4.1. The Simplified Problem We first consider this problem for the single- index 
sequence c^. 

To apply the Lagrange-Biirmann Inversion Theorem, we need to write the 
algebraic equation (|20p in the form G = z(j){G) for some formal power series (p. 
By simple arithmetic, one arrives at G = zcfyiG) with 

(note that 0(0) ^ 0). Then by Theorem mH 



Cd 



Since 



and 



(i-")'=E (!!)(-<' (30) 



n=0 ^ 



«~i:r''''W. m 



^ ' n=0 ^ 
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then Cauchy multiplication gives 



'»"-'"a^)'^|("'r>"G-i-,.)<-)^-'- 



._(^_^)d-iy^in + 2d-iy. d!(-l)"2" 



-^ n\{2d-l)l {d - I - n)l{n + ly. 



T1=0 

d-1 



(n + 2d-l)! d!(-l)" 2 



^ ^ ^^ (2d-l)! (d-l-n)!(2)„n! ^ ^ 

since (2)„ = 2 • 3 • • • (2 + n - 1) = (n + 1)!, and (-1)"(1 - d)n = Jl~%, gives 

=d(-^ \'i-^ V (^^ + 2rf-l)!(l-rf)«2" 
^ ) 2^ (2d-l)! (2)„n! 

(2d)„(l-d)„2" 



d-l 



=rf(-i)'^^^E 



n=0 
\d-l 



(2)„ n! 

=rf(-l)'^-i2Fi([l-rf,2d];[2];2). (33) 

We utilize (|29|) and Euler's formula: 

2Fi([a, 6]; [c]; z) = (1 - z)-'^2Fi (\a,c- b]; [c]; -^"j (34) 

(found in Gardner, ct. al. |10j . for instance) for the Gaussian hypcrgcometric 
function 

2Fi([a,6];[c];z) = f;^A^ (35) 

where (x)„ = x{x + 1) • • • (x + 7i — 1) is the Pochhammer symbol. With z = 2, 
c = 2, a = 1 — d, and 6 = 2 — 2d, (P^ and ([M]) give the hypergeometric function 
stated in the following: 

Theorem 4.2. The number of combinatorial invariants for S^ is 

Cd^2Fi{[2-2d,l-d];[2];2). (36) 

This series has a finite number of terms due to 1 — d e Z^ being non-positive for 
d>l. 

4.2. The Complete Problem We first find [z'']G(z,t). Following the same method 
as before, we note 

l + (l-(^ + ^^))G 

^^^^ = ii-{t+t^)Gr ' ^^^^ 

and by Theorem 14.11 we again have 
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Again utilizing Cauchy multiplication, we have 

n=0 \ / \ / 

and rewriting gives 

n=0 



^^ (2)„r7,! 



Recognizing that this is the right-hand side of ([M]). with a=l — d, 6 = 2 — 2(i, 
c = 2, and z ^ t + P, this simplifies beautifully to 

Cd{t)=2Fi{[l-d,2{l-d)];[2];t + t^). (41) 

By expanding (t + t^)", we get 

Theorem 4.3. The number of combinatorial invariants of dimension q for S^^ is 

c., = [t1.fi([l-d,2(l-d)];[2];i + t^) = |:i^^^Mi^f " ). (42) 

^ (2)„n! V? - V 

5. Asymptotic Behavior 

The next characteristic one might want to examine is the asymptotic behavior of 

the sequences c^ and Cd,q- 

5.1. The Simplified Problem For c^, we wish to utilize the relation 

1 



limsup^lcdl, (43) 

where R is the radius of convergence of the associated generating function. We first 
check that i? > by use of the Implicit Function Theorem. We have 



P(G, z) = G^ - G' + 


:^)- 


" 4 " 


= 


(44) 


dG (0,0) 4 








(45) 



so it follows that R > 0. We next determine R exactly in order to make a more 
precise asymptotic estimation for Cd- The polynomial P{G) behaves nicely; it has 
three roots, counted with multiplicity. That is, G has three branches, of which we 
are interested in the one going through (z, G) = (0, 0). The radius of convergence 
of the power series expansion of G about z = will be determined by the first 
singularity on the branch of interest. In this light, we wish to find all of the places 
the branches intersect, i.e. where the discriminant is zero. Calculation gives 
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so the discriminant A(P(G)) = when z = and z = ii±|^. When z == 0, P(G) 
has a double root at G = 1/2, and a single root at G = 0. The branch we arc 
interested in is the one where G = 0. So G is analytic in a neighborhood of z = 0, 
and the radius of convergence will be determined by the z with smallest modulus 
where there is a singularity. There are two numbers to check; we start with the 
one with smallest modulus, "'^"'^"2 ■ Using Sturm's Theorem, it can be shown that 
P{G) has three distinct real roots in the interval ]0, "l|^ [. At z = ii^|^, P{G) 
has a single root at G = ^( — 1 + \/5) and a double root at G = ^(3 — \/5). Since 
|(3 — VS) < ^(— 1 + -v/S), we can conclude that the branch we are interested in 
intersects the point {z,G) = C^Z^"^ , i(3 - v^)). Hence, R = ii^f^. Now since 
Cd is positive and increasing, equation (|43p becomes 



R 
and by the definition of a limit, we can conclude 



— = lim ^, (47) 



<'/«>' ^(ir^)'^'"^"'' 



cj ~ (1/R)« = — = » (1L09)^ (48) 



5.2. r/ie Complete Problem We consider the discrete probability distribution 

TTd., = ^. (49) 

Cd 

The goal is to show that this probability distribution is asymptotically of some 
known form, for example normal, when d — > c». We wish to use the same method 
as in [7]; compare with their Theorem 5. 

Theorem 5.1. The distribution iTd.q has mean /i^ and variance a^ that satisfy 

fid '^ nd, a1 ^ \d, (50) 

where k and A are algebraic numbers. The laws in each case are asymptotically 
normal for d — > 00 . 

Proof. The method of the proof is taken from [7] and also draws on the references 
[9], [3], and [8]. 

Consider the characteristic function of TTd{q), normalized by shifting the mean 
by fid and the variance by cr^: 

fd{t) = e-'^'^^l'^^Gd U^I^A , (51) 

where Gd(t) = ^^^^ is the probability generating function. The goal is to show that 
the normalized characteristic functions /d(i) converge pointwise to the characteristic 
function e~* 1"^ of the standard normal. Then the limit law follows by Levy's 
Continuity Theorem. 
We look at 

P{,z, i, G) = (< + t^fG"^ - 2{t + t^)G^ + {l-z + z{t + t^))G - z. (52) 
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By previous discussions, Go = (3 — \/5)/4 is a double root for P{p,l,G), where 
p = ii-5v5 ^ gQ i^y Weierstrass preparation, there is an analytic factorization 

P{z,t,G) = {{G - Gof + mi{z,t){G - Go) + m2{z,t)) ■ H{z,t,G), (53) 

where H{z,t,G) is analytic near (p, l,Go), il(p, l,Go) 7^ 0, and ttii and 7712 are 
analytic at (p, 1). Then applying the quadratic formula, 

G{z,t) - Go = i (-mi(z,i) ± ^mi{z,ty - 47n2(z,i)) ■ (54) 

We determine our branch of interest (as we did in Subsection 13. ip by the following. 
Consider first {z,t) restricted to < z < p{t) and < t < 1. Since G{z,t) is 
real there (G has real coefhcients), the discriminant mi{z,t)^ — 47712(2, t) must 
also be real and non-negative. Furthermore, since G is increasing in z for fixed t 
(coefficients of G are non- negative) , and since the discriminant vanishes at p{t) and 
is hence decreasing in z, then we need to take the minus sign, i.e. 

G(z,i) - Go = i (-mi{z,t) - v/77ii(z,t)2 -47^2(2,0) ■ (55) 

Set D{z,t) := 7Jii{z,t)'^ — Ara2{z,t). Consider the resultant 

R{z, t) = Resuhc {P{z, i, G), ^P(z, t, G)\ , (56) 

which is a polynomial whose restriction i?(z, 1) has, by the discussion in Section 
15.11 a simple root at p = ii-5V5 ^ gy |-j^g Implicit Function Theorem, this root lifts 
to an analytic branch p[t) of an algebraic function, for t near 1: 

i?(p(t),t)=0, p{l)^p. (57) 

Therefore, the function D{z, 1) has a simple real zero at z = p, and we can factorize 

D{z,t) = {p{t)-z)K{z,t), (58) 

for some analytic K satisfying K{p,l) ^ Q. So the uniform family of singular 
expansions for G{z,t) near {p{t),t) takes the form 



Giz,t) = co{t) + ci(i)v/l - z/pit) + 0(1 - z/p{t)), (59) 

uniformly with respect to t for t in a small neighborhood of 1, and with p(i), co(t), 
and ci{i) analytic at 1. If G{z,t) satisfies equation (|59p . we let z = z/p{t) and 
apply the results in |8] to obtain 

cd{t) := [z'']G{z,t) = (p(i))^''[z'^]{cl(i)^/T^I+0(l -z)} 

ci(^)(p(0)~'^ ,-, , ^^,-1 ,^,, .„r,^ 

The probability generating function is therefore a so-called "quasi-power," i.e. 
Gd{t) satisfies 
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If ([61]) holds, then for d — > oo and close enough to w = 1, 

G.«)~(^)'. (62) 

so for fixed t, letting d — 7> oo, 

^(0) 



fd{t) ^ e-^'^^*/"^ 



r{it/ad) 



(63) 



where r{s) = p{e'^)- Set y{s) = log(r(s)/r(0)), expand y(s) in a Taylor series, and 



use 



'■(0) '^ 

r(s) 



exp(dlog(r(0)/r(s))) to get 

fd{t) ^ exp{-zAidi/f7d - d {ity'{0)/ad - t^y"{0)/2aj + O(tV^d))}- (64) 

By calculation on (|5T|) . 

A*, = G',(l) ~ -dp'(l)/p(l) = -dy'(O) and 



^i = G'.'(i) + ..(1 - ..) ^ -d '"%''^'^ -{m))^ -''"^'^- ^''^ 



We obtain 

Ut) ^ e-*'/2 (66) 

for alH. □ 

By differentiation of (|5T|) and with help from MAPLE's algeqtodiffeq function in 
the GFUN package [12] . we can calculate 



so 



-11 + 5^5 87-39V5 702p(l) + 716p-(l) 
P(l) = ^ > P (1) = 2 ' P^^)"" Zqo ' (^^-^ 

n=^ = !^±±^ « 1.145898036, (68) 

p(l) -2 



and 



4^ - 4^ + f 4^V - ''"^''^ - 0.484597133. (69) 

p(i) p(i) Vp(i)y 10 ^ ' 



6. Alternative Counting Problems 

There are several variations on this enumeration problem that one might consider 

interesting. We discuss a few of these now. 

6.1. Distinguishing between Real and Complex Dimension It has been suggested 
that one might want to count the complex dimensions (corresponding to the 
number of analytic invariants in H for a vector field ^p) and the real dimensions 
(corresponding to the number of analytic invariants in R+ for a vector field ^p) 
separately (see [2] and [6] for definitions of the analytic invariants). This would 
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correspond in the complete problem to a triple-indexed sequence Cd^s,h, where s 
is the number of complex analytic invariants and h is the number of real analytic 
invariants. We follow the same method as in Section IX^ Let 

{"s"s„_i.s,/i if string does not start with parentheses 

"(s"s2a,si,/ii"s)"s6,s2Ji2 Si+ S2^ S, hi+h2 + l = h 
"[■5"s2a.si./ii"s]"Sb,S2./i2 Si + S2 + 1 = «, hi + 112 = h 

(70) 
be our string generating algorithm for 2a + b+2 = n. We define a recursion equation 
for Pn,s,h in the usual way, and arrive at the generating function G{z,t,'w) = 
'Yli'^d,s,hz'^'w^t^ satisfying (unsurprisingly) 

{t + wfG^ - 2(t + w)G'^ + {l-z + z{t + w))G - z = 0. (71) 

One can use methods as before to determine properties of C(i^s,h- 

6.2. Enumeration of Moduli Classes The combinatorial classes we have focused 
on distinguish between vector fields which are conformally conjugate by a rotation 
by a {d — 1)"' root of unity. Some may consider the problem of enumeration 
of classes in moduli space (i.e. without the labelling of the scparatriccs) more 
valuable. In order to study this, we utilize Polya theory for the cyclic group 
isomorphic to the group generated by e^^^/id~i)_ y^^, g^j-g interested in the sequence 
Cd = [z'^~^]C{G{z)) where G{G{z)) is generating function for the number of 
combinatorial classes for degree d vector fields moduli rotation by a (d — 1)'^* root 
of unity. Polya theory gives 

C(G(.)).i:Mlog^-l^. ,72) 

where </>(fc) is the Euler totient function. This relies on prime factorization, so it 
does not seem likely that one can find a closed-form expression for Cd as a function 
of d as of date. 
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